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Topological phase transitions between chiral and helical spin tex- 
tures in a lattice with spin-orbit coupling and a magnetic field 
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Abstract -We consider the combined effects of large spin-orbit couplings and a perpendicular 
magnetic field in a 2D honeycomb fermionic lattice. This system provides an elegant setup to 
generate versatile spin textures propagating along the edge of a sample. The spin-orbit coupling 
is shown to induce topological phase transitions between a helical quantum spin Hall phase and a 
chiral spin-imbalanced quantum Hall state. Besides, we find that the spin orientation of a single 
topological edge state can be tuned by a Rashba spin-orbit coupling, opening an interesting route 
towards quantum spin manipulation. We discuss the possible realization of our results using cold 
atoms trapped in optical lattices, where large synthetic magnetic fields and spin-orbit couplings 
can be engineered and finely tuned. In particular, this system would lead to the observation of a 
time-reversal-symmetry-broken quantum spin Hall phase. 



j Introduction. — The discovery of the first topological 
state of matter, the integer quantum Hall (QH) effect [1], 
has revolutionized our understanding of condensed mat- 
ter: a two-dimensional (2D) electron gas subjected to a 
perpendicular magnetic field is gapped and insulating in 
the bulk, but features chiral edge states that carry quan- 
tized charge currents. The quantization of the Hall con- 
ductance gh — v(e 2 /h) is directly related to the number 
of these chiral edge states z/, which is a sum of topological 
Chern numbers [2]. The proposal by Kane and Mele [3], 
that a new topological state of matter could be realized 
in graphene by taking into account spin-orbit interactions, 
has brought about the quantum spin Hall (QSH) effect, in 
which the edge currents carry spin but no charge. This 
spin transport is realized by helical edge states, namely 
spins traveling in opposite directions along the edge [cf. 
Fig. 1 (a)], and is protected by a Z2 topological index 
[3]. Contrarily to the integer QH effect, which is rooted 
in a time-reversal symmetry (TRS) breaking perturbation 
(i.e., a magnetic field), the QSH effect is a consequence 
of the so-called intrinsic spin-orbit coupling (ISO). The 
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latter acts as opposite "fluxes" [4] on the two spin com- 
ponents. Thus, the system effectively exhibits two QH 
phases, which get transformed into each other by time re- 
versal, giving rise to helical edge states. 

The two topological states of matter described above 
produce either charge or spin edge currents that are 
dissipation-free and robust against disorder, due to their 
topological character. In conventional semiconductors, 
which exhibit the integer QH effect, the spin texture of 
the edge states is fixed by the competition between the 
cyclotron energy and the Zeeman spin splitting produced 
by the external magnetic field. In particular, the observa- 
tion of Hall plateaus at high odd fillings v requires spin- 
resolved Landau levels which can be achieved by enhanc- 
ing the Zeeman energy, e.g., by tilting the sample [5] or 
by engineering materials with large Lande g factors [6] . In 
the limit of extremely large Zeeman energy, one eventually 
expects to observe spin-filtered QH states, where all the 
current-carrying edge states have the same spin orienta- 
tion [cf. Fig. 1 (b)]. 

An interesting question is whether topological phase 
transitions between helical QSH and chiral QH phases, 
exhibiting distinct spin textures on the edge, could be 
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driven in a system combining the effects of both a mag- 
netic field and spin-orbit couplings. In this Letter, we ad- 
dress this fundamental problem by solving the Kane-Mele 
model [3] in the presence of an external magnetic field. Al- 
though it is an interesting model from a theoretical point 
of view, its physical realization as a model for graphene has 
strong limitations: the topological phases emanating from 
the Kane-Mele model cannot be observed experimentally 
with graphene because of its extremely small ISO coupling 
[7]. The physics of this model can however be studied 
within an alternative and promising framework offered by 
cold atoms trapped in optical lattices. Indeed, control- 
lable Abelian and non-Abelian synthetic gauge fields can 
mimic the effects of magnetic fields and SO couplings in 
a highly versatile experimental framework [8-13]. In ad- 
dition, a Zeeman splitting term can also be easily engi- 
neered [14]. These recently developed technologies offer 
the unique possibility to explore the regimes of extremely 
large magnetic field and spin-orbit couplings. Moreover, 
cold-atom systems allow us to consider quantum phases 
in the absence of disorder, in which case the QSH state 
should remain stable even in the presence of TRS break- 
ing perturbations. 

Until now, most of the proposals for realizing QH lat- 
tice models with cold atoms are based on the spinless Hal- 
dane or Hofstadter models [15,16]. In these configurations, 
one expects to generate spin-filtered QH states, since only 
one spin component (i.e., internal atomic state) is con- 
sidered in the dynamics. On the other hand, an SU(2) 
generalization of the Hofstadter model, using more inter- 
nal states, has been introduced in Ref. [11] to engineer 
the QSH phase. However, despite their ability to either 
exhibit the QH or QSH phases, the lattice models of Refs. 
[11,15,16] do not fully capture the interplay between quan- 
tum Hall physics and the lifting of the spin degeneracy, 
present in condensed-matter systems due to Zeeman and 
spin-orbit effects. 

In this Letter, we show that a two-component lattice 
combining Zeeman and spin-orbit terms constitutes an 
ideal playground to produce versatile topological insulat- 
ing states with rich spin textures. It is shown that an 
optical-lattice realization of this system could display re- 
markable phase transitions between chiral and helical spin 
structures, not envisaged so far. The synthetic magnetic 
field and spin-orbit terms are generated by inducing a 
gauge field in the atomic cloud that influences all the hop- 
ping terms [8,10-12,15], whereas the Zeeman splitting fil- 
ters the spins. Different combinations of these terms will 
provide us with novel topological phase transitions that 
modify the spin texture characterizing the edge states. In- 
terestingly, these transitions are driven by the spin-orbit 
couplings, while the Fermi energy and external magnetic 
field are held fixed. Furthermore, we show that this setup 
allows to manipulate the spin orientation of an isolated 
edge state, thus generating richer spin structures than in 
standard topological insulating phases. 
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Fig. 1: Spin-textures and quantum Hall phases: Spin- 
up and spin-down components are represented by blue and red 
balls, respectively, (a) The QSH phase: the helical edge states 
are formed by counter-propagating spins of opposite sign, (b) 
The spin- filtered quantum Hall phase: the chiral edge states 
are constituted of identical spins with same orientation, (c) 
The spin-imbalanced QH phase: the chiral edge states feature 
several spins of different nature such that ^ N±. (d) The 
chiral Rashba phase: a single chiral edge state is formed with 
arbitrary spin direction, represented in purple. 

Results. — We consider the Kane-Mele model, which 
describes spin- 1/2 fermions in a honeycomb lattice, sub- 
jected to a uniform perpendicular magnetic field B = 
Bl z . The presence of the magnetic field modifies the 
graphene model of Ref. [3] in two fundamentally differ- 
ent manners. First, the magnetic field acts on the charges 
through the minimal coupling p — eA, where A is the 
gauge potential, and enters the model through the so- 
called Peierls phases [17]. Then, an additional coupling 
leads to a constant Zeeman splitting g\±BB<j zl where g 
is the Lande factor, \±b is the Bohr magneton and a z is 
the Pauli matrix acting on the spins. The tight-binding 
Hamiltonian of interest for our work then takes the follow- 
ing form 

U = -tJ2 ^ lk c{ci + iA SO Yl »k,ie [ * lk c{a zCl 

(k,l) «/c,Z» 

+ iA R ^ lk c{(o- x d M ) z ci + 27r$A z ^ c\cr z c k . (1) 

(k,l) k 

Here, = are the annihilation (creation) 

operators on the lattice site I for fermions with spin com- 
ponents up and down. The first two terms in Eq. (1) 
correspond, respectively to nearest (NN) and next-nearest 
(NNN) neighbour hopping on the honeycomb lattice, with 
the associated amplitudes t and Aso The second term 
is due to the ISO coupling and is characterized by spin- 
dependent amplitudes, with v^l = ±1 according to the 
orientation of the path connecting the NNN sites k and I. 
Here, we set H = e = 1 and we use the hopping parame- 
ter t and the bond length a = 1 as the energy and length 
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E* n = Mn+^n+l+^-l^n-l, 

v _ ( 2A S o^ sin(27r$(n + \) + k) + 27r$A z <7 ;2 tl - iA R 6- y 

tl + iARO^ -2Aso^ sin(27r^>(n + ~) + fc) + 27r^>A Z (3-; 
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, >/3. 1 A 

7± = ^-y^ + 2 a v ( 2 ) 



units, respectively. The Peierls phases, (j>ik = i Jj A • dZ, 
act on both the NN and NNN hopping terms and are gen- 
erally expressed in terms of the dimensionless parameter 

= SqB/2tt, i.e., the number of magnetic flux quanta 
per plaquette, where Sq = 3V3/2 is the area of a pla- 
quette. The third term corresponds to a Rashba SO cou- 
pling, which acts as a spin-mixing perturbation between 
NN sites: it tends to misalign the spin by rotating it to- 
wards the plane of the sample [3]. In this term, Ar is 
the coupling strength, & is the vector of Pauli matrices 
and dki is the vector between sites k and I. The last term 
in Eq. (1) is the constant Zeeman splitting, where we set 
Az = SI^bSq 1 . The potential realization of the model (1) 
with ultracold atoms trapped in optical lattices allows us 
to envisage the case of extremely large magnetic fluxes 

~ 0.1 — 1, where the spectral gaps A ~ 0.1 — It lead 
to Hofstadter fractal structures [17, 18]. In contrast to 
condensed-matter experiments (in which Az is basically 
fixed by the Lande factor), optical lattice setups allow us 
to tune every parameter individually to arbitrary values. 
In this framework, exotic configurations could be realized, 
e.g., featuring large magnetic flux and spin-orbit couplings 
while keeping the Zeeman splitting low. Furthermore, the 
non-interacting limit can be easily reached through Fesh- 
bach resonances [19]. 

In the following, we discuss novel topological phase tran- 
sitions and the manipulation of isolated edge states, which 
are obtained by tuning the Zeeman splitting and the SO 
terms. Although our calculations are performed for the 
honeycomb lattice, our predictions do not rely on this spe- 
cific geometry: Any system exhibiting a Dirac-like disper- 
sion relation, large Zeeman splitting and spin-orbit inter- 
actions will exhibit a similar behaviour. 

The spectrum and topological phases stemming from 
Hamiltonian (1) have been widely studied for two limiting 
cases. In the absence of magnetic flux, <I> = 0, a finite ISO 
coupling opens a QSH gap at half-filling, with magnitude 
A = 6v/3Aso [3]. In the spinless limit Az = Aso = Ar = 0, 
a finite magnetic flux leads to a multi-gap QH system 
[20]. In this Letter, we explore the combined effects of 
spin-orbit couplings and magnetic flux on these QH and 
QSH gaps. Since the honeycomb lattice has two sites per 
unit cell, A and £?, the single-particle wave function can be 
written as ^ n ,m = (^ m , ^n,™)? where (n, m) are integers 



labeling the unit cells [21], and where the spin components 
are implicit. Using the Landau gauge, we set Sfr nyrn = 
exp(i/cm)\I/ n and the single-particle Schrodinger equation 
takes the form of a generalized Harper equation, Eq. (2). 

Topological phases are distinguished by the quantum 
currents carried by gapless edge states, which are obtained 
by computing the spectrum E = E(k) on an abstract 
cylindrical geometry [22]. The transport coefficients are 
equally evaluated through the computation of topological 
invariants, the Chern numbers and Z2 index, which are 
associated with the bulk energy bands [2,3]. 

Let us first discuss the fate of the QSH phase as a uni- 
form magnetic field is progressively turned on. The main 
effect of this external field is to break TRS and there- 
fore to allow scattering processes between the counter- 
propagating edge states characterizing the QSH phase [3]. 
Therefore, although the bulk energy gap and the non- 
trivial topological index associated with the QSH phase 
at half-filling are preserved for a finite magnetic flux, the 
TRS violation destroys the topological phase in the pres- 
ence of disorder. Such phases, referred to as "weak QSH" 
(WQSH) phases or TRS-broken QSH phases, are thus triv- 
ial in real materials, but they would still be robust in op- 
tical lattices where disorder is absent. We note that the 
robustness of the WQSH has recently been investigated in 
Ref. [23] , for a system featuring a constant exchange term. 
Naturally, the magnetic field also creates additional spin- 
degenerate QH phases away from half- filling. For suffi- 
ciently large magnetic fields, the WQSH gap at half-filling 
is destroyed (e.g., when Aso = 0.05£ and Az = O.lt, this 
gap closing occurs at <I> « 0.14) and only the standard 
QH gaps remain. We highlight that the WQSH phase de- 
scribed here relies on the existence of a bulk energy gap, 
characterized by a non-trivial Z2 index, which guarantees 
the robustness of this phase for finite magnetic flux (in the 
absence of disorder). 

We now investigate an alternative and more interest- 
ing scenario, starting from spin-degenerate QH phases 
in the absence of Zeeman splitting and spin-orbit cou- 
pling. The energy spectrum for = 1/3 shows four 
spin-degenerate QH phases, with the Hall conductivities 
o~h = {2,— 2, 2, —2} in units of the conductivity quan- 
tum. Now, let us turn on the Zeeman splitting: by lifting 
the spin-degeneracy, the Zeeman term separates the QH 
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phases and eventually generates spin-filtered QH phases 
[cf. Fig. 1 (b) and Fig. 2 (b)]. When the Fermi energy is 
located inside such a spin- filtered QH bulk gap, a two- 
component Fermi gas will exhibit edge currents trans- 
ported by a single species only. Besides, the Zeeman term 
creates a WQSH phase at half- filling. It is remarkable 
that such a helical topological phase is produced in the 
absence of ISO coupling, which highlights the great simi- 
larity between the Zeeman and the ISO terms in Eq. (1) 
[24]. However, this Zeeman-induced QSH phase is weak 
because of the TRS violation [cf. discussion above] . Start- 
ing from this configuration, we progressively add the ISO 
coupling. As shown in Fig. 3 (a), the ISO opens new 
gaps in the energy spectrum, the topological character 
of which highly depends on the Zeeman splitting. For 
Az = 0.5£ and Aso = 0.35t, a WQSH gap opens at 
E w Lit [cf. Fig. 3 (a)]. Increasing the ISO term fur- 
ther leads to the closure of the gap, depicted in Fig. 3 
(b), and eventually to a topological phase transition. As 
represented in Fig. 3 (c), the gap is now characterized by 
a strong spin-imbalanced QH phase [cf. Fig. 1 (c)]. In- 
deed, the edge is populated by chiral edge states with 
different spin components: = 1 ^ N± = 2. In this 
configuration, when the Fermi energy E-p lies in the gap 
at E « 1.3t, the Hall conductivity and the spin-Hall con- 
ductivity differ, since gh = 3 and a s H = — 1 in units of the 
charge and spin conductivity quanta, respectively. Note 
that the chiral character of these current-carrying edge 
states guarantees the robustness of this quantum phase 
against external perturbations. In this sense, the ISO 
coupling has driven a transition from a weak QSH to a 
strong spin-imbalanced topological phase, while the Fermi 
energy and magnetic field were kept constant. Conse- 
quently, by adjusting the Zeeman splitting and the ISO 
coupling, one is able to tune the spin-imbalance ^ N± 
characterizing the current-carrying edge states. Recently, 
one-dimensional spin-imbalanced atomic gases have been 
under investigations [25], and lead to the realization of 
exotic pairing mechanisms. Our 2D model features spin- 
imbalanced states that are confined to the ID edge of the 
system, but also possess a chiral nature. We note that 
spin-imbalanced QH phases have already been observed 
in solid-state laboratories: In semiconductor heterostruc- 
tures without ISO coupling (e.g., GaAs/AlGaAs), a transi- 
tion from v = 2 (i.e., = N± = 1) to the spin-imbalanced 
phase v = 3 (i.e., = 2 and N± = 1) is driven by tuning 
the magnetic field. On the contrary, in our model the ISO 
interaction drives the transition, while the magnetic field is 
held constant. Furthermore, the variation of the magnetic 
field in semiconductors causes a phase transition between 
two chiral states, whereas in our study the ISO drives a 
transition between the WQSH phase, which is helical, to 
the chiral spin-imbalanced QH phase. Thus, the interplay 
between spin-imbalance, chirality and controllable inter- 
actions should lead to rich physics not envisaged so far. 

The Rashba SO coupling has dramatic consequences on 
the system described by Eq. (1), as it induces spin mixing. 



In the presence of a magnetic field, the Rashba term lifts 
the spin degeneracy of the original QH phases and even- 
tually opens new topologically trivial gaps. The Rashba 
term therefore destroys the topological nature of the QH 
gaps in certain regions of the spectrum (cf. Ref. [3] without 
magnetic field). In addition, the Rashba coupling partially 
polarizes the edge currents associated with the QH phases. 

A remarkable spin-manipulation process can be envis- 
aged by combining the effects of the Rashba spin-mixing 
perturbation and the Zeeman splitting. As already men- 
tioned above, the Zeeman splitting generates spin-filtered 
edge states in the QH gaps. In the example illustrated in 
Fig. 2 (b), a single spin-down particle propagates along the 
edge of the sample when E F « — 2.25t, for $ = 1/3 and 
Az = OM. As shown in Fig. 2 (c), one can then progres- 
sively increase the Rashba coupling and control the spin 
orientation (a^) = (^loyj^) of this specific and unique 
chiral edge state. This process indeed modifies the band 
structure, as well as the spin textures associated with the 
edge current, leading to the chiral Rashba phase illustrated 
in Fig. 1 (d): the spin components of this single edge state 
are now <j z = —0.5 and a± = 0.5, where a± is the spin 
component in the direction perpendicular to both the di- 
rection of propagation and the applied magnetic field. The 
canting of the edge-state spin is a hallmark of the compe- 
tition between the Rashba term, which tends to rotate the 
spin direction towards the xy plane, and the Zeeman term, 
which tends to align the spins in the vertical direction. For 
Ar = t, a large but realistic value for a Raman-induced 
tunneling amplitude, this single edge state can be reached 
by adjusting the Fermi energy to the value E? ~ — 3.25t, 
as shown in Fig. 2 (c). Interestingly, the spin direction 
of the edge states within this gap highly depends on the 
Fermi energy [26]. Thus, the spin direction of the edge 
states can be tuned by varying the filling factor. We em- 
phasize that this effect is notably different from the con- 
stant spin orientation that would result from the Zeeman 
effect in a tilted magnetic field. We finally stress that the 
variation of Ar does not close the topological bulk gap, and 
therefore, its associated Chern number remains unaltered. 
In this sense, there is no topological phase transition in- 
volved in this manipulation of the edge-state spin. 

Discussion and conclusions. — The experimental 
realization of the topological phase transitions presented 
here constitutes an appealing target for cold-atom labora- 
tories. In principle, one could realize the model described 
by Eq. (1) by trapping ultracold fermions in a honeycomb 
optical lattice. Indeed, this geometry has been recently 
realized experimentally, but we stress that simpler geome- 
tries with additional pseudo-spin degrees of freedom would 
lead to similar effects [12]. The main elements that have to 
be added to these lattice systems are the uniform magnetic 
field and spin-orbit couplings, which constitutes a difficult 
but appealing task. Both effects have already been engi- 
neered in cold gases confined to a single trap, by using Ra- 
man transitions between internal atomic states [8]. How- 
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0=1/3 +X z (=0.5t) +A 2 (=0.5f)+A R (=0 



Fig. 2: Generating a chiral Rashba phase on the edges: The energy spectrum E — E(k) is represented as a function of 
the momentum, in a cylindrical geometry with zigzag edges, (a) Several QH phases are generated by a magnetic flux $ = 1/3. 
When the Fermi energy lies in the gap at E « — £, a single spin-degenerate edge state propagates on each edge, (b) The 
inclusion of the Zeeman splitting lifts the spin degeneracy and creates spin-filtered QH phases at E — ±2.25£ and E — ±3£. 
At half-filling, the Zeeman splitting opens a "weak" QSH gap. (c) The addition of the Rashba term then allows for flipping 
the spin of the edge state to an arbitrary orientation. When the Fermi energy Ef « — 3.25£, the single edge state has spin 
components a z ~ —0.5 and a± « 0.5 (cf. orange dot). Here, the bulk energy bands are depicted in dark blue. 



3 Quantum spin Hall phase Q Gap closing point Q Spin-imbalanced QH phase 




Fig. 3: Phase transition from the QSH to the spin-imbalanced QH phase: The energy spectrum E — E(k) is 
represented as a function of the momentum, in a cylindrical geometry with bearded edges. Here we set $ = 1/3 and Az = 0.5t. 
Spin up [resp. down] states are represented by up [resp. down] triangles and the two colours designate the two edges of the 
cylinder where the states are localized, (a) The weak QSH phase is characterized by counter-propagating spins with opposite 
sign, (b) The gap closes for Aso ~ 0.4285£. (c) The spin-imbalanced QH phase: each edge features one spin-up and two 
spin-down states propagating in the same direction. Here, the bulk bands are depicted in black. 



ever, different methods are required for generating such 
fields in atomic lattices, as have already been suggested in 
Refs. [11,15,27]. In particular, we note that such a method 
has been recently described for the honeycomb lattice [28] . 
These proposals could be adjusted to generate the Rashba 
and ISO terms, as well as the Peierls phases. Moreover, 
the Zeeman splitting in Eq. (1) could be easily produced 
by static magnetic fields [14]. In this elaborate but fruitful 
configuration, the parameters Aso, Ar, and Az are finely 
tuned by external fields. This high control over the pa- 
rameters allows to drive the transitions proposed in this 
work. A difficulty encountered in cold-atom realizations 
of topological insulating phases is the presence of confin- 
ing traps, which will generally alter the bulk bands and 
destroy the edge-state structures [16]. To overcome this 
drawback, it is possible to either design sharp walls or to 
create interfaces in the system [11]. We also stress that the 



topological phases presented in Figs. 2 and 3 are protected 
by bulk energy gaps A ~ 0.1 — It, thus requiring achiev- 
able temperatures T ~ 10 nK for their observation. It is 
worth mentioning that the possibility to tune the inter- 
particle interactions with Feshbach resonances, combined 
with controllable Zeeman splitting and synthetic Rashba 
SO coupling, would lead to a cold-atom realization of a 
topological superconductor exhibiting Majorana fermions 
[29]. 

Finally, let us comment on the possibility to explore our 
results in a condensed-matter context, in which achiev- 
able magnetic fields typically belong to the low-flux regime 
(<£ <C 1). Recently, materials featuring strong ISO cou- 
pling and large Zeeman splitting, such as Hg(Cd)Te quan- 
tum wells, have been developed in order to detect the QSH 
[30,31] and spin-resolved QH states at low fields [6]. More- 
over, these setups allow experimentalists to study the com- 



p-5 



N. Goldman et al. 



petition between the Rashba SO and the Zeeman splitting 
[32]. In principle, the effects discussed in this Letter, which 
result from the combination of spin-orbit couplings and 
magnetic fields, could be obtained in these arrangements 
for realistic values of the model parameters. For exam- 
ple, manipulating the spin orientation of spin- filtered QH 
edge states could be achieved by tuning the electric field 
responsible for the Rashba coupling and by varying the 
Fermi energy within the bulk gap (cf. Fig. 2). However, 
these systems are unsuitable for detecting the topological 
phase transitions between helical WQSH and chiral QH 
phases. Indeed, a finite magnetic field breaks TRS and 
should therefore radically destabilize the WQSH phase in 
the presence of disorder: although the bulk gap and its 
non-trivial Z2 index might be preserved, the magnetic 
field triggers scattering processes between the counter- 
propagating (helical) edge states. Therefore, the WQSH 
phase induced by the Zeeman splitting (cf. Fig. 2) should 
reduce to a normal insulator in a solid-state experiment 
[3]. Moreover, we note that the ISO interaction, which 
drives the WQSH-QH transition (cf. Fig. 3), cannot be 
varied in real materials. In conclusion, cold atoms sub- 
jected to synthetic gauge fields and Zeeman splitting will 
certainly provide a richer framework to study these in- 
teresting topological phase transitions and spin manipula- 
tion. 
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